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1.0 INTRODUCTION 

The condensation of vapor on the inside surface of a tube is an important 
heat transfer mechanism in many thermal devices. An example of such a device 
is the heat pipe. This is a liquid-vapor two-phase system, often in the form 
of a closed tube, which transfers heat by evaporation and condensation at 
opposing ends of the tube. The heat pipe has been used to advantage in many 
aerospace and -satellite systems, and is now being recognized as a potentially 
attractive heat transfer device for terrestrial applications. The intentional 
or accidental introduction of noncondensable gases into heat pipes signifi- 
cantly changes the condensation heat transfer characteristics of the heat pipe 
condenser. This unique behavior has been widely used for temperature control 
in gas-loaded heat pipes. This report describes experiments and analyses done 
to further establish the characteristics of these gas zones in the limit as 
the gas becomes increasingly soluble in the working fluid. 

Early analyses of the gas-controlled heat pipe assumed a very simple model 
.for position and effect of noncondensable gas within the condenser; that is, 
the noncondensable gas was assumed to be in the form of a plug with a sharply- 
defined vapor-gas interface. In time it was found that the flat- front model 
did not accurately represent. many systems. It was observed that axial heat 
transfer via the condenser wall caused the sharp front to become more diffuse 
and that the boundary was also elongated due to vapor diffusion into the stag- 
nant gas zone. The culmination of one-dimensional models including wall con- 
duction was the numerical program GASPIPE^^^ which allowed the user flexibility 
in modeling condenser heat transfer effects and yielded better agreement be- 
tween design calculations and experimental observation. 


\ 
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In this study a more exact one- dimensional condensation heat transfer 
model for insoluble gases has been developed and compared with experimental 
data. Modifications to this model to accommodate soluble gas behavior have 
also been accomplished, and the effects on gas front behavior demonstrated. 

In Section 2, Analysis, analytical models for condensation heat transfer 
are documented. In Section 3, Experimental Apparatus and Technique, a novel 
optical method used for measuring gas concentration profiles is outlined. 
Section 3 also presents experimental data. Section 4 is an interpretation of 
these data and Section 5 is the summary and conclusion. 
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2.0 ANALYSIS 

In this section models are developed for condensation heat transfer and 
dissolved gas transport in a gas-loaded heat pipe. The influence on condensa- 
tion of an insoluble noncondensable gas is developed in Section 2.1, while 
various effects of a soluble gas are explored in Section 2.2. The intention 
of these sections is to construct an analytical framework that can be used for 
comparison with experimental heat transfer data presented in Section 3.8, 
Experimental Results. 

2. 1 One-dimensional Model for Condensation in the Presence of an Insoluble 

Noncondensable Gas 

To determine the effect of gas solubility on condensation heat transfer 
in a heat pipe, it is worthwhile to develop an analytical model for condensa- 
tion in the presence of an insoluble gas so that a bench mark or reference 
model is available for comparison with data. An additional incentive for 
developing such a model evolved during execution of this study— the soluble 
gas models, in addition to being complex, were numerically unstable. These 
problems are discussed in more detail in Section 2.2, Soluble Gas Modeling. 

The coupled' heat transfer and mass transfer processes in the condenser 
have been described by two one-dimensional conservation equations accounting 
for heat and mass transfer in the vapor phase and heat transfer along and 
through the heat pipe wall. Both axial and radial pressure gradients are 
neglected, and composition and temperature are assumed uniform across the 
entire vapor cross section. The boundary conditions used are similar to those 
of Marcus,^^^ but certain simplifying assumptions have not been made. Specif- 
ically, the first and second spatial derivatives for vapor and wall tempera- 
ture are not assumed equal, allowing accurate representation of high thermal 
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resistance systems such as the glass heat pipe used in experimental phases of 

this study. The modeling described in this section is similar to modeling by 
( 2 ) 

Saaski / ' but various improvements have been incorporated which increase 
solution accuracy and reduce computational time. 


2.1.1 Vapor Phase Mass Transfer 

The intent of this modeling is to define heat transfer in the stagnant 
gas zone shown in Figure 2.1 and existing between the axial positions 
z< z<z„. A cylindrical heat pipe cross section is assumed, but the 
equations to be developed are not geometry-dependent. However, uniform wall 
and vapor cross-sectional areas are necessary. An immiscible gas/fluid combi- 
nation is also assumed. If axial and radial pressure gradients and thermal 
conduction are negligible in the vapor phase, and the heat pipe wall is uni- 
form in temperature along any circumferential line, then vapor migration into 
the gas zone can be modeled as diffusive flow into a stagnant gas layer. Con- 
servation of working fluid vapor (species a) at any cross section results in 
the coupled heat and mass transfer equation 


dz 


\^fg 


(2-1) 


where denotes molar flux and T^ and T^ represent the vapor and wall tempera- 
tures respectively. is the molar heat of vaporization and the vapor cross 
sectional area is A^. This expression is based on the vapor phase being sepa- 
rated from the wall by a liquid film of uniform conductance per unit length, 
G^. This constant thickness film is attributable to a uniform thickness capil 
lary wick upon the wall. In the case of reflux operation without a wick, the 
film thickness is a slowly varying function of axial position, but can often 



FIGURE 2.1. REPRESENTATIVE THERMAL PROFILES FOR A GAS-LOADED HEAT PIPE 


be approximated by a constant value. The surface of this film which is in 
contact with the vapor phase is at the vapor phase temperature; any interfacial 
resistances are assumed small. For an insoluble gas, there is no net gas move- 
ment and (})jj = 0 (b denotes the noncondensable gas). The molar flux of species 
a is given as 


< 1 ) 


a 


CD dX /dz 



( 2 - 2 ) 


so that 
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The ultimate goal is to obtain a differential equation which has as 
the primary dependent variable. With this objective in mind, the mole frac- 
tion Xg in a constant pressure system is expressible as 


and 



Ps = f(T,); 


dT. 


= F’; 


d^P. 


s _ 


F" 


dT. 


(2-4) 


(2-5) 


where P^ is a function of the vapor temperature T^, and P^^ is the total heat 
pipe internal pressure. The concentration-diffusivity product, can be 

expressed as a function of T^ through the Chapman-Enskog formula for dilute 
gases 

‘^“ab ■ ^'ITT (2-6) 


where B is a very slowly varying function of temperature, and T^ is expressed 
in absolute temperature. For practical purposes, B can be evaluated at a mean 
temperature and assumed to be a constant. 

When the expressions above are substituted into equation (2-3), a second- 
order nonlinear differential equation in T^ results 



_WsnL_ 



(2-7) 


The numerical technique used to solve this equation is discussed in Section 

2.1.3. 
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2.1.2 Axial Wall Conduction 

The heat pipe wall is assumed uniform in cross section and to have a 
temperature dependence on the z-coordinate only. The wall is coupled to the 
vapor phase through a film or layer with per- unit- length thermal conductance 
and to the heat sink through a similar conductance, G^.. In Figure 2.2, 
this series of thermal conductances is shown as both a nodal network and as 
related to the condenser-heat sink assembly. A steady-state heat balance on 
an infinitesimal wall slice leads to the differential equation used in the 
experimental phase of this study: 

where T^ is the heat sink temperature and A^ and are the wall cross- 
sectional area and thermal conductivity. 

2.1.3 Numerical Solution Method - One- dimensional Problem 

In the last two sections, second-order differential equations for T^ and 
T^ have been derived. Under steady-state operating conditions, both equations 
must be satisfied simultaneously at all axial positions. 

By formation of a nodal network as shown in Figure 2.2, it has been pos- 
sible to solve the coupled equations (2-7) and (2-8) by converting the equa- 
tions to finite difference form and using successive over-relaxation (SOR) to 
iteratively approach consistent profiles for T^ and T^. However, the approach 
has required long computational times and alternate methods for solving the 
equations were investigated. All failed but one. Most notable among the 
failures was the so-called "shooting technique" in which initial guesses for 
the temperatures and their derivatives at z = are adjusted until the 
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SYSTEM COMPONENT DESIGNATION 



VAPOR TEMPERATURE 
FLUID FILM CONDUCTANCE 
WALL TEMPERATURE 

INTERFACIAL CONDUCTANCE 
HEATSINK 


Gf 

Tv. 

Gj 

Ts 


EQUIVALENT FINITE DIFFERENCE NODAL NETWORK 



CROSS SECTION OF 
HEAT REJECTION ASSEMBLY 


FIGURE 2.2. FINITE DIFFERENCE MODEL FOR GAS-GAP CONDENSER/HEAT SINK ASSEMBLY 
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solutions at z = z^, obtained by integrating the differential equations across 
the zone, satisfy a zero derivative boundary condition at z^^. However, the 
solution for would invariably either rise or fall with great speed after a 
number of apparently stable increments in z. The nonlinear equation (2-7) has 
many of the mannerisms of a "stiff" equation. 

Successive over-relaxation is commonly employed as a point-iterative 
technique; that is, one point at a time is adjusted in value. A block relaxa- 
tion technique was 'defined which was stable and improved convergence rates. 
Taken individually, the finite difference forms of (2-7) and (2-8) are tri- 
diagonal matrices for T^^. and T^^-, if temperatures on the other nodal line are 
fixed. That is, in the nodal network of Figure 2.2, if the wall temperatures 
are held fixed at one iterative step, then the system of linear equations de- 
fining the axial vapor temperature profile is tridiagonal. A tridiagonal' 

( 3 ) 

matrix can be solved exactly using the Thomas algorithm.^ ' This exact set of 
vapor temperatures is then held constant, and a new line of wall temperatures 
is calculated using the Thomas algorithm. The new wall temperatures are held 
constant and a vapor profile is calculated again as in the first step, complet 
ing the iterative cycle. All theoretical profiles were generated- using' this 
block iterative technique. 

Prior to generating difference equations, the temperatures and axial co- 
ordinate system were nondimensionalized as follows: 


T-T_ 


m s 


G.z‘ 

1 


A K 
w w 


(2-9) 


where T^^ is the adiabatic vapor temperature. The vapor pressure f(T^) was 
found by fitting experimental vapor pressure data to the following form: 
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anf = 9 q + ^ (2-10) 

The vapor temperature T^ is in absolute degrees. The finite difference form 
of equation (2-7) can be configured into a tridiagonal matrix with each row 
given by 


+ + + 
ot. T - T + 6. T . + Y- T .,T = d. 

1 v» 1-1 *^1 v,i 'i v,i+l 1 


( 2 - 11 ) 


For intermediate points excluding the boundaries, the constants are defined 
by the following set of equations: 


a, = 1-E, 

(2-12) 

B,. = -(2 + F,) 

(2-13) 

Y, = lYE. 

(2-14) 

('f T„.i) v,1+l l.i-l) [l-x/ U 2 ] 

(2-15) 


(2-16) 


(2-17) 


(2-18) 


(2-19) 
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The finite difference equations for the wall were similarly derived and 
will not be presented. Boundary conditions used in modeling were 



1.0 

0.995 

1/(1 + G./G^) 



+ + 
z < 


+ + 
z = z^ 


+ + 
z = z„ 
0 


0.75 


+ 


z 


z! + 8.5 
0 


(2-20) 

( 2 - 21 ) 

( 2 - 22 ) 

(2-23) 


Standard methods were used to incorporate these boundary conditions into 
the finite difference equation sets; for example, the derivative boundary 
conditions were expressed as central divided differences. 

The dimensionless vapor temperature was calculated over the range 

-f -f -f 

z^ < z s Zq + 8.5. At a distance of 8.5 dimensionless units from the gas 
front, the gas concentration is essentially constant and calculations over a 
longer zone are not necessary and reduce solution accuracy in the area around 
Zq for a fixed number of increments. The axial gas zone was divided into 114 
increments. Wall temperatures were calculated over a wider range extending 
into the nonblocked condenser zone. This is done because the heat pipe wall 
in the gas zone acts as an extended fin, draining heat away from the warmer 
areas and depressing wall temperature beyond the actual gas front. 

A listing of the computer program used for modeling purposes is given in 
Appendix C. The code is in BASIC language and was programmed on a PDP-11 
minicomputer. Representative profiles are tabulated in Appendix B, and a 
comparison with experimental data is given in Figure 3.10. Iteration of 
temperature profiles was terminated when the largest change in wall 
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temperature per iteration was less than 0.001 in dimensionless units. This 
typically required 40 iterations--at the end of these iterations, the vapor 
temperature was often very slightly less than the wall temperature for 
» Zq. With a larger number of iterations this physically unrealistic 
situation is corrected, but the difference in profile shapes is quite negli- 
gible, and was not considered worth the additional computer run-time. Forty 
iterations required about 5 to 6 minutes of CPU time on the PDP-11. 

2,2 Soluble Gas Modeling 

An initial attempt to model condensation in a tube containing soluble 

gas was done using a von Karman-Pohl hausen integral technique. By assuming 

characteristic equations for the concentration and velocity profiles in both 

vapor and liquid phases, it was possible to generate coupled differential 

equations describing the two-dimensional concentration and flow fields for 

incompressible flow. The approach used was similar to that of Somogyi and 
(4) 

Yen, ’ The analysis developed will not be presented here in detail because 
the authors were unsuccessful in obtaining stable numerical solutions, 

A shooting technique was used from the interface at z = z^. A stable 
solution set was obtained until the differential integration approached the 
fully stagnant gas plug zone. Under these conditions, the fluid film on the 
heat pipe wall is saturated with a large amount of noncondensable gas, and 
the gas-saturated liquid moves toward the evaporator into regions of 
lower vapor-phase gas concentration. This supersaturation creates a large 
gas flux out of the fluid surface and, in fact, the integral techniques indi- 
cated total vapor phase fluxes that were directed radially inward. When this 
occurred, the solution became unstable and was not correctable with smaller 
steps in z. Because of these large radially inward gas fluxes, it was not 
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possible to obtain a zero derivative for concentration and zero value for the 
axial and radial velocity components at z - z^^. The integral method was 
abandoned and the hybrid 1D/2D method described in the following section was 
developed. 

2.2.1 Hybrid Soluble Gas Model 

For tubes on the order of 1 cm diameter and vapor pressures of 1 

atmosphere (typical heat pipe conditions), vapor phase diffusivity can be 

much higher than liquid phase diffusivity. Vapor phase diffusivities range 
‘ 2 

typically from 0.1 to 1.0 cm /sec, while liquid phase diffusivities of 1 to 
-5 2 

10 X 10 cm /sec are typical for gases in common working fluids. 

Therefore, if radial concentration gradients are neglected in the vapor 
phase but are accounted for in the liquid phase, it is possible to extend the 
analysis of Section 2.1 to include soluble gas effects for a wide range of 
physically relevant heat pipe designs. 

If interfacial resistance is neglected at the vapor/liquid interface 
as in Section 2.1, then the concentration of gas at the liquid phase inter- 
face is in the equilibrium with the gas concentration in the vapor phase, and 

= a Cgy (at liquid/vapor interface) (2-24) 

where is the liquid phase gas concentration, is the vapor phase gas 
concentration, and a is the Ostwald coefficient. A mass balance on the 
working fluid vapor yields the same expression as equation (2-1). However, 
the noncondensable gas is no longer stationary, and the equation (2-2) for 
(f), must be modified to 

a 
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'as 


1 - X. 


(2-25) 


and the conservation equation for coupled heat transfer and vapor transfer is 


A 

dx 


(♦as) 


^(Tv-Tv,) 


(2-26) 


To obtain the flux <|)jj in (2-25), transport of noncondensable gas within 
the surface liquid layer must be detailed. In the experiments described in 
Section 3, the heat pipe was wickless and vertically oriented with the 
evaporator at the bottom. The surface layer therefore is in the classic 
falling film category as investigated by Nusselt and others. If it is 
assumed that the dissolved gas has no effect on physical properties of the 
liquid and that diffusion of dissolved gas in the vertical direction is negli- 
gible compared to convective transport, then conservation of gas in the film 
is governed by 

V(y,z) ^ = D (2-27) 


where V(y,z) is the film velocity and is the diffusion coefficient of the 
gas in the liquid. The coordinate y is measured radially inward with respect 
to the heat pipe inner wall surface. 

Under conditions of laminar flow, negligible shear at the liquid-vapor 
interface, and negligible vapor density compared to liquid density, the 
velocity profile is 



(2-28) 
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where 


+ 

y 


= y/6 


(2-29) 


The fluid density and viscosity are p,^ and respectively. The fluid 
film thickness is 6. At the wall, V{0) is zero, while at the vapor/liquid 

I 

interface, V(l) is a maximum. 

Boundary conditions on Cg are 


o 

CQ 

II 

aC 

gv 

+ 1 
y = 1 

(2-30) 

dC „ 

gs. _ 
dy 

0 

+ 

11 

o 

(2-31) 


aC 

gvm 

Z = 0 y"^ 1 

(2-32) 


At z = z^, the dissolved gas concentration has reached a maximum value in 
equilibrium with the fully stagnant gas concentration, 

The liquid film thickness 6 changes with axial position because of 
condensation. By using the velocity distribution (2-28) and assuming a 
linear temperature profile across the film, film growth is obtained by 
integrating over the heat transferred from z = z, to z = z , where z z 
z^. The film thickness 6 is given by 


t 5 P 

w .15 

* 4 K 









(2-33) 


where t is the heat pipe wail thickness. 

For coupling to the finite difference equations it was necessary to 
include the wall thermal resistance between the fluid film and the node at 
the center of the wall. For metal -walled heat pipes this radial wall 
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temperature drop can be neglected, but it is of first-order importance in a 
glass heat pipe because of a much lower thermal conductivity. 

For numerical solution of the conservation equation (2-27), a dimension- 
less gas concentration and radial measure were defined. 


C 


+ 




(2-34) 


The gradually-changing film width 6 must be accounted for in the con- 
servation equation 


V(y^,z) • 


3z " 6 


dz 






(2-35) 


The conservation equation (2-35) with boundary conditions 2-30 through 
2-32 was solved by implicit finite difference methods similar to those 
described in Section 2.1. The gas phase concentrations used in boundary 
conditions 2-30 were obtained from a finite difference solution of the vapor 
phase equation (2-26). This solution was also obtained in a manner entirely 
analogous to the solution method in Section 2.1 except for inclusion of the 
new term in i .e. , 

a 

new term = (2-36) 


The axially-varying flux <j)j^ was obtained from the solution to the liquid- 
phase diffusion problem. At each axial position, a mass balance including 
both the vapor and liquid phases shows that the flux in the vapor phase 
must be balanced by a gas flux of opposite sign within the liquid phase. 
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The liquid-phase gas flux is obtained by numerically integrating the 

concentration profile. 


Vbi 


= ctTrD^6 




: /■' 

gvm / 

Jq 


V(y*) 


. + 
c dy 


(2-37) 


where is the heat pipe inner diameter. 

The coupled problem was .iteratively solved by alternately holding the 
solution fixed in one phase and solving the finite difference equation in the 
other phase. At the completion of each cycle, <J>|3(z) would be updated. When 
the concentration profiles in each phase remained constant with iteration, 
the coupled problem was considered solved. 



18 


3.0 EXPERIMENTAL APPARATUS AND TECHNIQUE 

3. 1 Alternate Methods for Measuring Vapor Phase Gas Concentration 

To compare experimental data with analytical models for condensation in 
the presence of a soluble gas, an approach is desired which indicates 
accurately both vapor phase composition and spatial variations in composition. 
Many previous investigators have inferred internal vapor phase dynamics from 
variations in heat pipe wall temperature, or from probe temperatures in the 
vapor phase. The first method has the obvious deficiency of not being a 
direct measurement and, in fact, it is easy to show that measured axial 
temperature variations are primarily related to axial wall conduction for 
most metal -walled heat pipes. It is very difficult to infer from external 
measurements an accurate picture of vapor phase mass transfer in the transi- 
tion zone between the freely condensing and fully stagnant zones. 

It is also difficult to accurately measure the condition of the vapor 
phase directly. Any probe inserted into the vapor generally has such a large 
thermal mass that profiles can be severely distorted. Thermocouple techniques 
or thermistors are generally difficult to integrate into an operating heat 
pipe, not only affecting the accuracy of temperature measurement through their 
own mass, but also fundamentally changing the flow distribution. 

If a small diameter thermocouple well is placed along the axis of a 
cylindrical heat pipe, boundary conditions for the Navier-Stokes equations 
are modified as: 

BEFORE PROBE: atr = 0, U = U, |^=0 

o’ 9r 



AFTER PROBE: 


at r = 0, U = 0, 
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where U represents the axially-directed velocity component. This is a very 
fundamental change in velocity distribution and, therefore, the measurement 
technique has significantly influenced the system being measured. 

3.2 Vapor Composition by Optical Spectroscopy 

As discussed, spatial variations in noncondensable gas concentration are 
difficult to measure within an operating heat pipe without obtaining a dis- 
torted picture of vapor composition. This difficulty has been overcome to a 
great extent by using optical spectroscopy. 

The basic concept is shown in Figure 3.1. A working fluid/gas combina- 
tion is selected in which the noncondensable gas is either strongly absorbent 
or optically transparent at a wavelength where the working fluid vapor is of 
opposite characteristic. Using a monochromatic light source, such as a laser 
or spectrophotometer, the outgoing beam from the source is shaped by a narrow 
slit perpendicular to the heat pipe axis. The heat pipe enyeJLope is glass, 
and light transmitted through the heat pipe impinges upon a sensitive photo- 
cell coupled to a picoammeter. The heat pipe/condenser assembly is mounted 
on a vertical drive so that the heat pipe can be scanned by the beam. 

Absorption of light by a gas is given by the familiar Lambert-Beer Law 

I = EXP(-BV,^) (3-1) 

where is the incident light intensity, p is an absorption coefficient, and 
is the distance of light travel through the absorbing medium. is the 
effective, geometrically-weighted average concentration along the path at any 
given axial position. At constant wavelength a photo-diode has a linear 
response to light intensity, so that I and can also be regarded as photo- 
detector current. Therefore, if p and are known, concentration of the 
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FIGURE 3.1. BASIC TECHNIQUE USED FOR MEASURING CONCENTRATION PROFILES 
OF OPTICALLY-ABSORBING GASES IN A REFLUXING TWO-PHASE SYSTEM 
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absorbing species, Cj^, can be found. 



The constant a now is a cell constant and assumes uniform geometry along the 
heat pipe length. 

This technique allows very accurate measurement of gas concentration 
versus axial position with minimum disturbance to the heat pipe vapor. 

Spatial variation along the heat pipe axis can be defined to any degree of 
accuracy by decreasing slit width, although practical limits arise in obtain- 
ing adequate light intensity to drive the photodetector. 

The technique also has some limitations. It is obvious that the 
technique is principally useful in defining the average concentration as a 
function of axial position; the beam in effect integrates any radial varia- 
tions. It is also not advisable to use any wick upon the heat pipe inner 
diameter, as it will optically scatter light. Therefore, this technique is 
most useful for studying condensation in the presence of noncondensable gases 
upon smooth vertical surfaces. Such surfaces are of considerable importance 
in two-phase heat transfer equipment; for example, many terrestrial appli- 
cations use heat pipes in vertical reflux, in exact agreement with the method 
optimum for testing. In addition, both vapor phase and condensate phase 
conditions are very similar to conditions in zero-g, where a wall wick and 
surface tension pumping replace the falling condensate film defined by viscous 
shear and gravity, 

3.3 Apparatus Description 


Table 3.1 summarizes the experimental system. 
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TABLE 3.1. EXPERIMENTAL TEST VEHICLE 
CHARACTERISTICS 


Heat Sink 
Type 

Length 

Diameter 

Optical slot width 
Optical slot height 
Heat pipe centering method 

Conductance of centering 
screws 

Heat Pipe 

Envelope type 
Wall thickness 
Nominal length 
Evaporator length 
Adiabatic length 
Working fluid 
Gas additive 

Optical Assembly 

Operating wavelength 

Effective beam size at 
heat pipe 

Nominal indicator current 
with no absorbing gas 

Background with light 
source shielded 

Resol uti on 


Gas-gap chiller block with optical 
slot; Aluminum 

15.24 cm 

1.4618 ±0.0013 cm (I.D.) 

0.318 cm 
13.20 cm 

Four 6-32 nylon screws at each end 
of sink 

0.012 W/K for each set of four 
screws 


1.395 cm O.D. Pyrex glass tube 
0.12 cm 

50.0 cm 

16.0 cm 
18.8 cm 
CCl^ 

Chlorine 


3550 A 

0.10 X 0.318 cm 

4.03 (10"^^) ampere 

<0.02 (10“^°) ampere 

±0.025 (10~^^) ampere on 10"^*^ scale 
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Figures 3.2a and 3.2b present perspective views of the glass heat pipe 
and aluminum heat sink. The heat sink, which was in the form of a split-shell, 
had an inner diameter of 1.416 ±0.0013 cm, and was 15.24 cm long. Four cool- 
ing channels of 0.63 cm diameter maintained the condenser isothermal. The 
optical transmission slits were 0.318 cm wide and 13.2 cm long. After the 
heat sink was machined, it was given a dark brown chromate conversion coating 
having a high absorptivity in the blue portion of the spectrum. The outside 
was then painted flat black. In Figure 3.2b, the light beam is visible in 
the optical slit. 

Before operation, the glass heat pipe of 1.345 cm O.D. was inserted into 
the inner hole, and centrally located using a feeler gauge and 8 radially 
mounted nylon 6-32 screws; four screws are at the top of the condenser and 
4 at the bottom, spaced at 90-degree intervals. These screws served to hold 
the heat pipe in the condenser with a minimum conductive loss, and establish 
a uniform gas-gap circumferentially. 

The heat pipe was 50 cm in overall length. The evaporator was 16 cm 
long. To maintain the evaporator and adiabatic sections as thermally isolated 
as possible, an aluminum walled Dewar, continuously pumped with a nitrogen 
trap and fore pump, surrounded these sections, and was sealed to the heat pipe 
at the condenser end with a Buna-N 0-ring. Electrical feed-throughs supplied 
current to the evaporator resistance heater through the opposite end of the 
Dewar. On the external surface of the Dewar, a uniformly-space resistance 
heater maintained the exterior of the Dewar within 1°C of the heat pipe adia- 
batic temperature. This was necessary to reduce thermal losses via the 0-ring 
contact at the condenser end, and via a mounting plate at the opposite end 
which centrally located the heat pipe in the Dewar. These elaborate pre- 
cautions were necessary since the heat pipe was operated as much as 75°C above 
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FIGURE 3.2a. PERSPECTIVE VIEW OF TEST SYSTEM USED IIM MEASUREMENTS OF AXIAL 
NONCONDENSABLE GAS CONCENTRATION DISTRIBUTIONS. BY USING AN 
OPTICALLY ABSORBING NONCONDENSABLE GAS, A PRECISE BEAM OF 
MONOCHROMATIC LIGHT CAN DEFINE CONCENTRATION PROFILES. 
VARIATIONS IN TRANSMITTED LIGHT INTENSITY ARE DIRECTLY RE- 
LATED TO CHANGES IN GAS CONCENTRATION 



FIGURE 3 2b FRONTAL VIEW OF CONDENSER AND HEAT SINK SECTIONS 
TRANSMITTED LIGHT BEAM IS VISIBLE NEAR THE BOTTOM 
OF THE HEAT SINK GROOVE. 
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laboratory ambient. The condenser was capable of moving vertically via a 

> 

friction drive mechanism so that any portion of the transmission slit could 
be intercepted by the detector beam. 

3.4 Gas and Fluid Selection 

Based upon a high absorption in the blue spectrum, chlorine was selected 
as the noncondensable gas. A typical absorption spectrum for Cl^ gas at 
1 atmosphere and 20°C is shown in Figure 3.3. An operating wavelength of 

o 

3550 A was selected for all optical tests. This wavelength is very close to 
the absorption peak of Figure 3.3. As a working fluid, carbon tetrachloride, 
CCl^, was chosen on the basis of compatibility with chlorine, no significant 
blue spectrum absorption for thin fluid films, and a high affinity for taking 
chlorine into liquid solution. The latter characteristic is essential to 
define the effect of solubility on gas zone behavior. The solubility of Cl 2 
in CCl^ and R-11 was experimentally measured using an optical absorption 
technique. Results of these experiments are shown in Figure 3.4. 

3. 5 Calibration: Optical System 

It was experimentally established that the logarithm of optical absorp- 
tion was not linearly dependent upon gas concentration. Experimental 
absorption measurements with known concentrations of chlorine in a 1.4 cm 
diameter cell established the following relationship between gas content in 
standard cm^ per cc and £n (Iq/I) at A = 3550 A and 20.8°C. 

Cg = 0.39815 An + 0.066075 (An(I^/I))^ An(^)<1.35 (3-3) 


PRECEDING PAGE BLANK NOT PEiMEQ 
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WAVELENGTH (MILLIMICRONS) 

FIGURE 3.3. THE ABSORPTION SPECTRUM OF 1 ATMOSPHERE OF CHLORINE GAS AT 
20°C AS A FUNCTION OF WAVELENGTH. AN OPTICAL PATH LENGTH OF 
APPROXIMATELY 1 CM WAS USED. MONOCHROMATOR BAND-WIDTH 
WAS NOT DETERMINED. 
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FIGURE 3.4. THE SOLUBILITY OF CHLORINE IN CCI^ AND R-11 (CFCI3) 

Cg = 0.29267 - 0.30541 S,r\ (1^/1) + 0.49546 (jln (1^/1))^*^ 

/^o\ (3-4) 

1.35 

The experimental data basis for these expressions is shown in Figure 3.5. 

At a higher temperature T, the concentration is 

Cg(T) = 1 - 1^0.0005441 £n{I^/I) (T-20.8)J • Cg(20.8°C) (3-5) 

The light source employed at all times was a 500-watt slide projector 
operated from a constant- voltage transformer. The wavelength selected for 

o 

operation, 3550 A, was produced with an EG&G monochromator. Model 585-22. 
Beam thickness at the heat pipe was chosen at 0.1 cm as a compromise between 
spatial resolution and photodetector resolution. The radiometer detector 




DENSITY OF CU GAS (std cc/cm'*) 
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head and indicator units were also manufactured by EG&G, and were Models 
580- 22A and 580- llA, respectively. 

3.6 Thermal Measurements 

Heat sink and heat pipe temperatures were measured using Type K thermo- 
couples. Adiabatic temperature was measured with a thermocouple centered in 
the adiabatic section on the heat pipe exterior wall. All temperatures were 
measured with Q.1°C resolution using a Doric DS-350 digital thermocouple 
indicator. 

3.7 Data Reduction 

At completion of a test series, data included heat pipe vapor core 
temperatures and heat sink temperatures, heat input data, and photometer 
current readings as a function of axial position in the heat pipe. 

Heat losses or gains through the Dewar surrounding the evaporator and 
adiabatic sections were negligible. However, the nylon locator screws 'at 
the condenser entrance form a direct heat leak to the heat sink, and this 
heat loss was subtracted from the gross heat input to obtain net heat de- 
livered to the condenser. Each set of nylon screws had a combined conduc- 
tance of 0.012 W/°K. 

Gas concentrations were calculated using Equations (3-3) through (3-5), 
with the sink temperature used as the characteristic temperature in equation 
(3-5). The reference photocurrent was based on the peak current immedi- 
ately in front of the gas front. For most data runs, this current was about 
4.03 X 10 ampere. Total gas content within the gas zone was calculated 
by numerically integrating the concentration profile with Simpson's Rule. 

Gas concentration beyond the farthest measureable condenser position (about 

13.7 cm) was assumed equal to the concentration at 13.7 cm. 
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3.8 Experimental Results 

Concentration profiles from experimental sequences using CI 2 /CCI 4 at 
sink temperatures of 14°C and 40°C are shown in Figures 3.6 and 3.7, 
respectively. The axial position scale is with reference to the condenser 
end nearest the evaporator; e.g., the blind end of the condenser is dis- 
placed 16 'em axially from the "front end" of the condenser. The continuous 
curves, which represent theoretical predictions, are discussed in Section 4. 
The raw and converted data are both given in complete tabular form in 
Appendix A. 

Vapor temperature response to heat input is displayed in Figure 3.8 
for both sink temperatures. Figure 3.9 presents the gas zone molar charge 
based on integration of the experimental concentration profiles in Figures 
3.6 and 3.7. The molar gas load is shown versus stagnant zone gas concen- 
tration to determine if gas is going into solution in the stagnant zone. 
Figure 3.10 presents two models for predicting noncondensafale concen- 
tration distributions compared with data for the 2.2 W high temperature 
CCl^ run. 

Analysis of the data presented here is discussed in the following 
section. Data Interpretation. 



IMONCOIMDENSABLE CONCEiMTRATION (std cm^ per cm' 



AXIAL CONDENSER POSITION (cm) 

FIGURE 3.6, A COMPARISON OF EXPERIMENTAL CONCENTRATION PROFILES WITH THE THEORETICAL ONE- 
DIMENSIONAL MODEL ASSUMING NO SOLUBILITY EFFECTS. GAS CONCENTRATION PROFILES 
WERE CALCULATED FROM OPTICAL ABSORPTION DATA, SINK TEMPERATURE = 14,°C 



NONCONDENSABLE CONCENTRATION {std cm per cm 



FIGURE 3 7 A COMPARISON OF EXPERIMENTAL CONCENTRATION PROFILES WITH THE THEORETICAL ONE- 
DIMENSIONAL MODEL ASSUMING NO SOLUBILITY EFFECTS. GAS CONCENTRATION PROFILES WERE 
CALCULATED FROM OPTICAL ABSORPTION DATA SINK TEMPERATURE = 40.°C 




VAPOR-SINK TEMPERATURE DIFFERENCE (°K) 



0 2 4 6 8 10 12 14 16 18 20 

HEAT TRANSFER RATE (WATTS) 


FIGURE 3.8. A COMPARISON OF CALCULATED TEMPERATURE-HEAT TRANSFER CHARACTERISTICS WITH EXPERIMENTAL 
DATA CALCULATED VALUES ARE BASED ON THE KNOWN VAPOR-SINK TEMPERATURE DIFFERENCE AND 
THE NUMERICALLY INTEGRATED TOTAL GAS CONTENT. THIS INFORMATION IS INPUT TO THE qNE- 
DIMENSIONAL CONDENSATION HEAT TRANSFER MODEL DESCRIBED IN ANALYSIS SECTIONS, AND A GAS 
PROFILE AND HEAT TRANSFER RATE ARE THEN CALCULATED. 
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STAGNANT ZONE GAS CONCENTRATION (std cm^/cm^) 


FIGURE 3 9 TOTAL AMOUWT OF NONCONDENSABLE GAS CONTAINED WITHIN STAGNANT GAS 
ZONE AS A FUNCTION OF STAGNANT ZONE GAS CONCENTRATION. GAS CONTENT 
OBTAINED BY NUMERICALLY INTEGRATING EXPERIMENTALLY-DETERMINED 
CONCENTRATION PROFILES. NOTE THAT THE VAPOR PHASE APPEARS TO LOSE 
CHLORINE GAS AS THE GAS CONCENTRATION INCREASES. 




DENSITY OF Cl, GAS (std cc/cm 



POSITION RELATIVE TO GAS ZONE LEADING EDGE (cm) 


FIGURE 3 10 TWO MODELS FOR PREDICTING NONCONDENSABLE CONCENTRATION DISTRIBUTIONS COMPARED WITH DATA 
FOR THE 2 2W HIGH TEMPERATURE CCt. RUN (SINK=40°C) BY MODELING VAPOR DIFFUSION AT THE LEAD- 
ING EDGE OF THE STAGNANT GAS ZONE, MUCH BETTER AGREEMENT WITH EXPERIMENTAL DATA RESULTS 
NO SOLUBILITY EFFECTS ARE APPARENT 




37 


4.0 DATA INTERPRETATION 

4.1 Comparison of Experimental Results With Theory 

The full complement of data profiles was simulated with the one- 
dimensional insoluble gas model and the results are shown in Figures 3.6 and 
3.7 as the solid curves. The experimental runs were simulated on the basis 
of the following data. 

1. The vapor-sink driving temperature difference 

2. The numerically-integrated experimental total gas content 

3. Physical property data and conductances. 

With these input data and the governing differential equations, the gas 
concentration profiles were calculated and the gas front positions determined. 
The analytical. and experimental data are in good general agreement, but the 
model rather consistently predicts a somewhat shorter gas zone. 

The general underestimation of gas zone length by the analytical model 
is not attributable to the governing differential equations but rather to a 
discrepancy between the optical absorption-derived gas density in the stagnant 
zone and the value calculated from the ideal gas law. This difference is 
summarized in Table 4.1, a comparison of heat transfer and gas concentration 
data with theory. Although the difference in concentration amounts on the 
average to only about it is the cause of the gas front shift. The heat 
transfer rate as a function of driving temperature difference, given in 
Figure 3.8 and Table 4.1, is also in agreement to within ±5%; again, the 
differences are primarily attributable to the discrepancy in gas zone length. 

In the low heat flux runs the optically-deri ved stagnant zone concen- 
trations are in better agreement with the ideal gas law, and the data can be 
compared with the model in greater detail, as in Figure 3.10. If axial vapor 
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TABLE 


Sink Tempera- 
ture (°C) 

14 nom. 


40 nom. 


.1. A COMPARISON OF EXPERIMENTALLY-DETERMINED VALUES 
FOR HEAT TRANSFER AND GAS CONCENTRATION WITH VALUES 
PREDICTED BY THE ONE-DIMENSIONAL ANALYTICAL MODEL 


Stagnant Zone 


Vapor-sink Gas Concentration 

Temperature Heat Transfer Rate (W) (std cm^/cm^) 


Difference (“O 

Expt 

Model 

t Dev 

Expt 

Model 

% Dev 

35.65 

2.09 

2.05 

-1.9 

.285 

.293 

2.8 

38.8 

4.58 

4.57 

-0.2 

.330 

.340 

3.0 

41.8 

6.95 

7.00 

0.7 

.377 

.388 

2.9 

45.05 

9.37 

9.54 

1.8 

.427 

.445 

4.2 

48.2 

11.95 

12.03 

0.7 

.482 

.508 

5.4 

51.45 

14.3 

14.56 

1.8 

.565 

.579 

2.5 

58.17 

18.8 

19.24 

2.3 

.691 

.740 

7.1 

23.23 

2.22 

2.01 

-9.5 

.330 

.323 

-2.1 

26.9 

4.67 

4.57 

-2.1 

.394 

.398 

1.0 

30.8 

7.13 

7.21 

1.1 

.474 

.485 

2.3 

34.75 

9.65 

9.72 

0.7 

.567 

.582 

2.6 

38.65 

11.92 

12.33 

3.4 

.660 

.690 

4.5 

42.85 

14.3 

14.9 

4.2 

.754 

.817 

8.4 

51.1 

18.8 

20.2 

lA 

>.993 

1.104 
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transport were completely neglected in the modeling, then the gas concen- 
tration would be defined by a fin- type axial temperature profile in the 
stagnant region. This profile has been superimposed on Figure 3.10 to show 
the influence of vapor diffusion into the stagnant gas. The concentration 
data of Figure 3.10 are correlated well by the one-dimensional insoluble gas 
model. Under most operating conditions, the gas concentration profiles did 
not show any gross effect of the highly soluble gas. 

There was, however, a definite "reduction" in noncondensable gas as the 
stagnant zone gas concentration increased, i.e., as the heat transfer rate 
was increased. It is very probable that this was due to gas being absorbed 
into liquid in the stagnant zone. At a sink temperature of 14°C the Ostwald 
coefficient for chlorine in carbon tetrachloride is 31.0, while it is 18.5 
at 40°C. Assuming no loss of gas from the gas zone, the total moles of gas, 
N^, must remain a constant. For a sharply defined gas zone in length, and 
for a trapped volume of liquid 

N„ = ( Y, + aV ) (4-1) 

But in an optical absorption experiment such as done here, only the gas-phase 
noncondensable is seen, and this observable amount is 

SvmVv = "o - 

Therefore, if stays constant, it can be predicted that the observed moles 
of gas will appear to linearly decrease with increasing stagnant gas concen- 
tration, and that the rate of decrease will be proportional to the Ostwald 
coefficient. Furthermore, in the limit of going to zero, all test runs 
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at different sink temperatures should intersect at = 0 and this will 
represent the total gas content, N^. 

Measured total contents shown in Figure 3.9 are in agreement with all 
predictions above. The linear rate at which gas is "lost" at 14°C is about 
1.9 times the rate at 40°C, and the Ostwald coefficient is about 1.7 times 
higher at 18°C than 40°C. In addition, intersection of the two linear gas 
profiles in Figure 3.9 occurs very close to zero gas concentration. There- 

3 

fore, the total molar gas content in the heat pipe is about 4.75 std cm . 

It is informative to now use (4-2) to calculate the amount of liquid 
contained as droplets and surface film within the stagnant zone. This calcu- 
lation indicates a liquid load of only 0.039 cm is sufficient to produce 
the virtual gas loss shown in Figure 3.9. 

This effect was the basis of another study' '' where it was shown that 
under certain conditions a liquid absorption matrix can be used as a gas 
reservoir for gas-controlled heat pipes. This technique was shown experi- 
mentally to reduce reservoir volume by a factor of more than five compared 
to standard gas reservoir designs. 

The high solubility of chlorine in CCl^ also produced another unique 
effect not normally encountered in gas-loaded heat pipes. In the tests at 
low sink temperature, there was an unusual increase in chlorine concentration 
in the gas front region. This is shown very well in the 18.8 watt run of 
Figure 3.6. Because the heat pipe is very nearly a constant-pressure system, 
this zone must actually be cooler than the heat sink. In the fully stagnant 
zone, the wall and liquid films are close to sink temperature; a higher con- 
centration of chlorine in an intermediate zone implies a local vapor pressure 
(and film temperature) that must be below sink temperature. Another related 
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and unusual manifestation of high power level operation was observed in 

similar experiments carried out with an R-11 heat pipe and chlorine gas. In 

this case, droplets continuously formed within the vapor phase and fell 

through the gas front zone until they impinged on upward-moving vapor. This 

created a stagnation zone of high droplet density. A time-lapse photograph 

of this phenomenon is shown in Figure 4.1. The vertical streaks are falling 

droplets, while the lower diffuse disk is the droplet stagnation zone. This 

photograph was taken through the optical slit of the test apparatus in Figure 

3.2a with strong back lighting. This phenomenon is occurring within a tube 
2 

of only about 1 cm cross-sectional area. 

At the 18.8 watt heat flux for each sink temperature, the gas front was 
also very unstable with large "temporal" fluctuations in the axial concen- 
tration profile. The position of the gas front shifted axially 0,15 to 
0.20 cm in an erratic fashion with a time constant on the order of 15 to 45 
seconds, and the vapor core temperature also fluctuated about 0.1°C. 

The 14°C profile is shown in greater detail in A of Figure 4.2, along 
with a theoretical profile based on the soluble-gas modeling developed in 
Section 2.2. The axial positions of the experimental and theoretical profiles 
have been adjusted so that the profiles overlap. The experimental data range 
is indicated by the bounds given as dotted lines in A of Figure 4.2, while 
the theoretical profile is displayed as a solid line. The effects of the 
high chlorine gas solubility on the previously discussed insoluble gas con- 
centration profiles were very modest, and the differences in profile would be 
indistinguishable in A of Figure 4.2, This was primarily because of the thin 
fluid films (on the order of 0.001 cm) associated with the vertical reflux 
mode.. 



FIGURE 4.1 A TIME-LAPSE PHOTOGRAPH OF DROPLET FORMATION SEEN 
AT 18 WATTS IN A GLASS HEAT PIPE CHARGED WITH CI2 AND 
R-11. THE DROPLETS ARE BEING PHOTOGRAPHED WITH 
STRONG BACKLIGHTING THROUGH THE OPTICAL SLOT 
SHOWN IN FIGURE 3 2a. THE DROPLETS FALL UNTIL RESIS- 
TANCE IS ENCOUNTERED FROM UPWARD-MOVING VAPOR, 
CREATING THE DIFFUSE DISC-SHAPED AREA BELOW THE 
FALLING DROPLET STREAK PATTERNS ONE-SECOND 
ELAPSED TIME. 
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FIGURE 4.2. CONCENTRATION PROFILES AND CALCULATED RADIAL MOLAR FLUX DENSITIES 
FOR THE 19 WATT, 14°C SINK EXPERIMENTAL RUN. NOTE BAND-STRUCTURE DUE 
TO FLUCTUATIONS IN GAS CONCENTRATION. 
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The unsteady gas front behavior is believed to be related to radial free 
convection. If a Grashof number is calculated based on the tube radius and 
the vapor-sink temperature difference, a value of about 170,000 is found. 

For free convection between parallel vertical plates the critical Grashof 
number is only 8000. If significant radial concentration gradients were 
present in the transitional zone from free condensation to stagnation, then 
the radial Grashof number could be high enough to initiate a convective 
motion, and this unstable convective motion would explain the concentration 
fluctuations in A of Figure 4.2. This could also explain the cold zone indi- 
cated by absorption measurements at the 14°C sink temperature. A possible 
scenario leading to this behavior will now be discussed. 

In B of Figure 4.2, the calculated radial fluxes of vapor (4>y|^) and non- 
condensable gas at the liquid-vapor interface are shown. The vapor 

condensation rate shows a maximum at a point 0.2 cm into the gas zone, and 
decreases monotonically beyond this point. Chlorine gas is simultaneously 
being evolved from the film between 0 < z^. < 0.35 cm, and dissolved into the 
film from 0.35 cm to the end of the transition zone. Both the vapor and 
chlorine radial flux densities are small but finite at 1.05 cm, the position 
of maximum chlorine concentration in A of Figure 4.2. If the vapor phase gas 
concentration rises and falls in response to the disordered mass transfer 
occurring in the first 1/2 cm of the zone, chlorine will be absorbed and 
desorbed from the liquid film. From Figure 3.4, the heat of solution of Clg 
in CCl^ is about 14,000 Joules/g-mole, as calculated from the slope of log a 
versus 1/T, while the heat of vaporization of CCl^ is about 31,000 Joules/ 
g-mole in this temperature range. Because the film has 20 to 30 times as 
much Cl 2 per unit volume as the vapor phase, it is possible that gas desorp- 
tion from the thin film under semistagnant conditions could lead to a net 
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cooling of this zone, and a relative maximum in chlorine vapor phase concen- 
tration. However, because the modeling and experimental phases of this study 
were basically one-dimensional approaches, it is not possible to substantiate 
the mechanisms described above, and these processes should be regarded as 
speculatory at this time. 

4.2 Model Validity 

The one-dimensional models developed in this study assume that axial 
transfer of energy is predominated by diffusive mass transfer and that the 
sensible heat capacity of the diffusing vapor is negligible compared to the 
latent heat capacity. These assumptions are validated in Table 4.2, which 
compares axial transport of heat via gas-phase conduction, ordinary diffusion 
and thermal diffusion. The comparison was based on the theoretical profile 

TABLE 4.2. THE RELATIVE IMPORTANCE OF VARIOUS ENERGY 
TRANSFER PROCESSES FOR A REPRESENTATIVE 
CONCENTRATION PROFILE 

Axial Heat Flux Density 
Mode (W/cm^) 


Conduction 

0.0073 

Thermal di f fusion^ 

0.0004 

Latent heat, ordinary diffusion 

0.092 

(2) 

Sensible heat, ordinary diffusion' ' 

0.015 


^Assuming a thermal diffusion ratio of 0.05. 

( 2 ) ' 

^ Based on the total vapor-sink temperature difference. 
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in A of Figure 4.2. The comparison was done at the axial position with maxi- 
mum vapor temperature rate-of-change versus z^. This position was = 0.38 
cm. Ordinary diffusion is shown in Table 4.2 to be the primary mode of energy 
transfer, with heat transfer by conduction accounting for less than 8% of 
the total thermal flux, and thermal diffusion accounting for less than 0.5% 
of the energy transfer. The initial assumptions pertaining to energy trans- 
port modes are therefore generally satisfied by model solutions. However, 
the energy associated with vapor phase sensible heat was neglected, yet it 
could account for up to 15% of the axial energy flux. Inclusion of sensible 
heat would be an appropriate basis for a more accurate model of gas front 
concentrati on prof i 1 es . 
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5.0 SUMMARY 

0ne-diiTiensiona1 models for condensation in the presence of insoluble and 
soluble gases have been developed for a tubular geometry. The models have 
been compared with experimental data derived from Clg/CCl^ vertical reflux 
heat pipes and are in good agreement with measured concentration profiles at 
low heat flux densities. At high heat flux densities, the gas fronts become 
unstable and exhibit a turbulent and fluctuating concentration profile. This 
turbulent behavior was associated with unusual effects, including an apparent 
depression of heat pipe wall temperature below sink temperature and the for- 
mation of liquid droplets in the vapor phase. 

The experimental concentration profiles were determined using a new 
technique having several advantages over conventional thermocouple or thermis- 
tor temperature measurement methods. By selecting a noncondensable gas which 
strongly absorbs light over discrete portions of the spectrum, the concentra- 
tion of noncondensable gas as a function of axial position was determined 
directly from optical absorption measurements in a glass heat pipe. Axial 
resolution is improved over temperature measuring methods and the vapor-phase 
flow field is undisturbed. 

Apart from the effects at high power level, gas solubility did not have 
a significant influence on gas concentration profiles. This was predicted 
by the soluble gas condensation model, and was physically attributable to the 
relatively thin fluid films (on the order of 0.001 cm) associated with the 
vertical reflux mode. The short time constants for liquid-phase diffusive 
equilibrium precluded any significant effect on gas front shape. This may 
not be true for systems with thicker fluid films or lower liquid phase dif- 
fusivities, and each case must be considered individually. 
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Liquid films and droplets within the stagnant portion of the gas-rich 
zone were found to create a virtual sink for noncondensable gas. That is, 
as vapor phase gas concentration was increased, the absorption of increasing 
amounts of gas within this liquid was detected as a 10 to 20% loss in total 
moles of noncondensable gas in the vapor phase. This phenomenon could be 
used in gas-controlled heat pipes to increase thermal control capability or 
reduce gas reservoir size by dissolving the noncondensable gas into a volu- 
metrlcally efficient liquid reservoir in lieu of the common gas-bulb 
reservoir. 
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APPENDIX A 

COMPREHENSIVE TABULATION OF OPTICAL 
ABSORPTION AND THERMAL DATA 



Run A1 
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Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


49.65°C 

19.2°C. 

2.087 W 
4.342 std cm^ 
4.03(x 10“^°) A 


Axial position (cm) 


Current (x 10 


-10 


)A 


Concentration 
(std cm^/cm^) 


1.294 

3.960 

.0070 

1.344 

3.940 

,0091 

1.394 

3.830 

.0205 

1.444 

3.730 

.0313 

1.494 

3.580 

.0482 

1.544 

3.420 

.0674 

1.594 

3.260 

.0878 

1.644 

3.100 

.1095 

1.694 

2.920 

.1358 

1.744 

2.810 

.1529 

1.794 

2.690 

.1727 

1.844 

2.550 

.1972 

1.894 

2.410 

.2235 

1.944 

2.310 

.2436 

1.994 

2.240 

.2583 

2.044 

2.185 

.2703 

2.094 

2.150 

.2781 

2.144 

2.130 

.2827 

2.194 

2.120 

.2850 

2.244 

2.110 

.2873 

2.294 

2.110 

.2873 

2.344 

2.090 

.2920 

2.394 

2.075 

.2955 

2.494 

2.080 

.2943 

2.594 

2.080 

.2943 

2.694 

2.090 

.2920 

3.194 

2.080 

.2943 

4.194 

2.080 

.2943 

5.194 

2.090 

.2920 

6.194 

2.090 

.2920 

7.194 

2.090 

.2920 

8.194 

2.120 

.2850 

9.194 

2.120 

.2850 

10.194 

2.120 

.2850 

11.194 

2.120 

.2850 

12.194 

2.130 

.2827 

13.194 

2.120 

.2850 
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Run A2 


Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


52.8°C 
14°C 
4.584 W 
4.299 std cm^ 
4.03(x A 


Axial position (cm) 


Current (x 10 ^*^)A~ 


Concentration 
(std cm^/cm^) 


1.694 

4.080 

-.0049 

2.194 

4.080 

-.0049 

2.694 

4.070 

-.0039 

2.944 

4.070 

-.0039 

3.094 

4.030 

.0000 

3.194 

4.020 

.0010 

3.294 

3.970 

.0060 

3.344 

3.890 

.0142 

3.394 

3.780 

.0259 

3.444 

3.630 

.0425 

3.494 

3.390 

.0711 

3.544 

3.170 

.0998 

3.594 

2.910 

.1373 

3.644 

2.740 

.1643 

3.694 

2.610 

.1865 

3.744 

2.510 

.2045 

3.794 

2.380 

.2294 

3.844 

2.300 ■ 

.2457 

3.894 

2.180 

.2714 

3.944 

2.080 

.2943 

3.994 

2.020 

.3088 

4.044 

1.970 

.3212 

4.094 

1.930 

.3315 

4.144 

1.925 

.3328 

4.194 

1.890 

.3420 

4.244 

1.890 

.3420 

4.294 

1.890 

.3420 

4.344 

1.890 

.3420 

4.394 

1.890 

.3420 

4.494 

1.890 

.3420 

4.594 

1.930 

.3315 

4.694 

1.930 

.3315 

5.194 

1.940 

.3289 

6.194 

1.930 

.3315 

7.194 

1.930 

.3315 

8.194 

1.940 

.3289 

9.194 

1.940 

.3289 

10.194 

1.940 

.3289 

11.194 

1.940 

.3289 

12.194 

1.940 

.3289 

13.194 

1.935 

.3302 



Run A3 
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Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


55,8°C 
14.05°C 
6.949 W 
4.22 std cm^ 
4.03 (x 10"^°) A 


Axial position (cm) 


-TO 

Current (x 10 )A 


Concentration 
(std cin^/cm^) 


1.694 

4.080 

-.0049 

2.194 

4.080 

-.0049 

2.694 

4.060 

-.0030 

3.194 

4.040 

-.0010 

3.694 

3.940 

.0091 

4.194 

4.020 

.0010 

4.694 

4.030 

.0000 

4.794 

4.030 

.0000 

4.894 

3.990 

.0040 

4.994 

3.910 

.0121 

5.044 

3.820 

.0216 

5.094 

3.630 

.0425 

5.144 

3.430 

.0662 

5.194 

3.120 

.1067 

5.244 

2.880 

.1419 

5.294 

2.680 

.1744 

5.344 

2.540 

.1990 

5.394 

2.380 

.2294 

5.444 

2.270 

.2519 

5.494 

2.170 

.2736 

5.544 

2.065 

.2979 

5.594 

1.940 

.3289 

5.644 

1.870 

.3474 

5.694 

1.800 

.3668 

5.744 

1.780 

.3725 

5.794 

1.740 

.3842 

5.844 

1.710 

.3932 

5.894 

1.730 

.3872 

5.944 

1.730 

.3872 

5.994 

1.730 

.3872 

6.094 

1.730 

.3872 

6.194 

1.730 

.3872 

6.694 

1.760 

.3783 

7.194 

1.780 

.3725 

8.194 

1.780 

.3725 

9.194 

1.790 

.3696 

10.194 

1.780 

.3725 

11.194 

1.780 

.3725 

12.194 

1.790 

.3696 

13.194 

1.764 

.3771 



Run A4 
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Vapor core temperature: 59°C 

Sink temperature: 13.95°C 

Heat input: 9.369 W 

3 

Calculated gas content: 4.142 std cm 

Zero-absorption current: 4.02(x 10"^^) A 


Axial position (cm) 


Current 


[x 10“^^)A 


I 


2.194 

4.080 

3.194 

4.020 

4.194 

3.970 

5.194 

3.980 

5,694 

3.990 

5.794 

3.940 

5.894 

3.990 

5,994 

3.990 

6.094 

3.980 

6.194 

3.970 

6.294 

3.940 

6.394 

3.900 

6.444 

3.830 

6.494 

3.770 

6.544 

3.640 

6.594 

3.370 

6.644 

3.130 

6.694 

2.820 

6.744 

2.610 

6.794 

2.410 

6.844 

2.260 

6.894 

2.120 

6.944 

2.030 

6.994 

1.890 

7.044 

1.835 

7.094 

1.730 

7.144 

1.650 

7.194 

1.590 

7.244 

1.570 

7.294 

1.560 

7.344 

1.560 

7.394 

1.560 

7.694 

1.580 

8.194 

1.600 

9.194 

1.610 

10.194 

1.610 

11.194 

1.610 

12.194 

1.610 

13,194 

1.600 


Concentrati on 


[std cm^/cm^) 


-.0059 
.0000 
.0050 
.0040 
.0030 
.0081 
.0030 , 
.0030 
.0040 
.0050 
.0081 
.0122 
.0195 
.0259 
.0403 
.0726 
.1042 
.1502 
.1853 
.2224 
.2529 
.2838 
.3051 
.3408 
.3557 
.3859 
.4105 
.4300 
.4367 
.4401 
.4401 
.4401 
.4334 
.4267 
.4234 
.4234 
.4234 
.4234 
.4267 
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Run A5 


Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 

62.3°C 
14.1°C 
11.952 W 
4.046 std cm^ 
4.04(x 10"^^) A 

Concentration 

Axial position (cm) 

Current (x 10~^^)A' 

(std cm^/cm^) 

2.194 

4.120 

-.0078 

3.194 

4.040 

.0000 

4.194 

3.980 

.0060 

5.194 

3.960 

- .0080 

6.194 

3.930 

.0111 

7.194 

3.980 

.0060 

7.494 

3.990 

.0050 

7.594 

3.980 

.0060 

7.694 

3.940 

.0101 

7.794 

3.780 

.0269 

7.844 

3.640 

.0424 

7.894 

3.370 

.0747 

7.944 

3.090 

.1120 

7.994 

2.760 

.1621 

8.044 

2.475 

.2122 

8.094 

2.275 

.2521 

8.144 

2.120 

.2862 

8.294 

1.980 

.3199 

8.244 

1.873 

.3478 

8.294 

1.780 

.3738 

8.344 

1.680 

.4037 

8.394 

1.590 

.4326 

8.444 

1.500 

.4637 

8.494 

1.425 

.4915 

8.544 

1.410 

.4972 

8.594 

1.376 

.5106 

8.644 

1.367 

.5142 

8.694 

1.370 

.5130 

8.794 

1.370 

.5130 

8.894 

1.400 

.5011 

8.994 

1.405 

.4992 

9.094 

1.406 

.4988 

9.194 

1.415 

.4953 

10.194 

1.480 

.4709 

11.194 

1.450 

.482d 

12.194 

1.450 

.4820 



Run A6 
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Vapor core temperature: 65.65°C 

Sink temperature: 14.2°C 

Heat input: 14.263 W 

•3 

Calculated gas content: 3.92 std cm 

Zero-absorption current: 3.98{x lO'^^) A 


Axial position (cm) 

2.194 

3.194 
■ 4.194 

5.194 

6.194 

7.194 

8.194 

8.394 

8.494 

8.594 

8.694 

8.794 

8.894 

8.944 

8.994 
9.044 
9.094 
9.144 

9.194 
9.244 
9.294 
9.344 

9.394 
9.444 

9.494 
9.544 

9.594 
9.644 

9.694 
9.744 

9.794 
9.844 

9.894 

9.944 

9.994 
10.094 

10. 194 

10.694 

11.194 

11.694 

12.194 

12.694 

13.194 

13.694 


Current (x lO'^^lA- 

3.830 

3.820 

3.840 

3.920 

3.940 

3.980 

3.970 
3.940 
3.880 
3.920 
3.930 
3.930 
3.850 
3.730 
3.580 
3.380 

2.980 
2.680 
2.360 
2.180 

1.970 
1.844 
1.740 
1.630 
1.550 
1.460 
1.400 
1.320 
1.270 
1.220 

1.210 

1.185 

1.210 

1.210 
1.215 
1.215 
1.230 
1.280 
1.305 
1.120 
1.300 
1.300 
1.276 
1.230 


Concentration 

(std cm^/cm^) 

.0154 

.0165 

.0144 

.0061 

.0040 

.0000 

.0010 

.0040 

.0102 

.0061 

.0051' 

.0051 

.0133 

.0262 

.0431 

.0671 

.1213 

.1687 

.2275 

.2654 

.3150 

.3481 

.3778 

.4117 

.4382 

.4702 

.4930 

.5253 

.5468 

.5694 

.5741 

.5859 

.5741 

.5741 

.5717 

.5717 

.5648 

.5424 

.5317 

.5295 

.5338 

.5338 

.5442 

.5648 
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Run A7 

Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 

72.27°C 
14. 1°C 
18.802 W 
3.534 std cm^ 
4.03(x 10"^°) A 

Concentrati on 

Axial position (cm) 

Current (x 10"^°)A 

(std cm^/cm^) 

2.194 

4.300 

-.0256 

3.194 

4.100 

-.0069 

4.194 

4.000 

.0030 

5.194 

3.900 

.0132 

6.194 

3.850 

.0184 

7.194 

3.820 

.0216 

8.194 

3.600 

.0459 

9.194 

4.020 

.0010 

10. 194 

3.980 

.0050 

10.294 

3.980 

.0050 

10.394 

4.030 

.0000 

10.494 

4.000 

.0030 

10.594 

3.880' 

.0152 

10.694 

3. 780' 

.0259 

10.794 

3.750' 

.0291 

10.894 

3.350 

.0761 

10.994 

2.950 

.1313 

11.094 

2.500 

.2064 

11.194 

1.890 

.3420 

11.294 

1.480 

.4696 

11.394 

1.300 

.5407 

11.494 

1.160 

.6052 

11.594 

1.090 

.6413 

11.694 

.965' 

.7121 

11.794 

.926 

.7371 

11.894 

.907 

.7499 

11.994 

.916 

.7438 

12.094 

.945 

.7247 

12.194 

.960 

.7152 

12.694 

1.000 

.6910 

13.194. 

1.000' 

.6910 

13.694 

1.000 

.6910 



Run 81 
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Vapor core temperature: 63.35°C 

Sink temperature: 40.12°C 

Heat input: 2.221 W 

3 

Calculated gas content: 4.498 std cm 

Zero-absorption current: 4.03(xl0”^^) A 

Concentrati on 

Axial position (cm) Current (x 10~^^)A (std cm^/cm^) 


1.294 

3.920 

.0111 

1.694 

4.000 

.0030 

1.794 

4.030 

.0000 

1.894 

4.030 

.0000 

1.994 

4.010 

.0020 

2.094 

4.020 

.0010 

2.194 

3.980 

.0050 

2.294 

3.890 

.0142 

2.344 

3.810 

.0227 

2.394 

3.660 

.0392 

2.444 

3.540 

.0530 

2.494 

3.380 

.0724 

2.544 

3.230 

.0918 

2.594 

3.030 

.1195 

2.644 

2.925 

.1350 

2.694 

2.820 

.1512 

2.744 

.2.705 

.1699 

2.794 

2.620 

.1844 

2.844 

2.520 

.2022 

2.894 

2.470 

.2115 

2.944 

2.420 

.2210 

2.994 

2.370 

.2307 

3.044 

2.320 

.2408 

3.094 

2.270 

.2511 

3.144 

2.220 

.2617 

3.194 

2.180 

.2704 

3.244 

2.165 

.2737 

3.294 

2.140 

.2793 

3.344 

2.120 

.2838 

3.394 

2.080 

.2930 

3.494 

2.035 

.3037 

3.594 

1.970 

.3196 

3.694 

1.980 

.3171 

4.194 

1.970 

.3196 

5.194 

1.970 

.3196 

6.194 

1.970 

.3196 

7.194 

1.970 

.3196 

8.194 

1.940 

.3272 

9.194 

1.960 

.3221 

10.194 

1.970 

.3196 

11.194 

1.970 

.3196 

12.194 

1.990 

.3146 

13.194 

1.940 

.3272 

13.694 

1.930 

.3297 



Run B2 
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Vapor core temperature: 57.1°C 

Sink temperature: 40.2°C 

Heat input: 4.667 W 

3 

Calculated gas content: 4.421 std cm 

-10 

Zero-absorption current: 4.03(xl0 ) A 


Axial position 


Current 


Concentration 


'std cm^/cm^) 


1.694 

4.030 

.0000 

2.194 

4.060 

-.0030 

2.694 

4.060 

-.0030 

3.194 

4.060 

-.0030 

3.694 

3.970 

.0060 

4.194 

4.020 

.0010 

4.294 

4.020 

.0010 

4.394 

4.010 

.0020 

4.494 

3.990 

.0040 

4.594 

3.980 

.0050 

4.694 

3.910 

.0121 

4.744 

3.820 

.0216 

4.794 

3.680 

.0368 

4.844 

3.530 

.0540 

4.894 

3.360 

.0747 

4.944 

3.180 

.0981 

4.944 

3.020 

.1204 

5.044 

2.740 

.1634 

5.094 

2.610 

.1853 

5.144 

2.500 

.2050 

5.194 

2.355 

.2326 

5.244 

2.265 

.2509 

5.294 

2.180 

.2690 

5.344 

2.120 

.2824 

5.394 

2.070 

.2939 

5.444 

2.025 

.3045 

5.494 

1.980 

.3155 

5.544 

1.940 

.3254 

5.594 

1.'930 

.3280 

5.644 

1.890 

.3383 

5.694 

1.880 

.3409 

5.744 

1.840 

.3516 

5.794 

1.820 

.3570 

5.894 

1.800 

.3626 

5.994 

1.790 

.3654 

6.094 

1.780 

.3682 

6.194 

1.753 

.3759 

7.194 

1.730 

.3825 

8.194 

1.720 

.3854 

9.194 

1.730 

.3825 

10.194 

1.710 

.3884 

11.194 

1.710 

.3884 

12.194 

1.730 

.3825 

13.194 

1.700 

.3914 

13.694 

1.690 

.3944 



Run B3 
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Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


70.95°C 
40.15°C 
7.13 W 

4.337 std cm^ 
4.03(xl0"^°) A 


Axial position (cm) 


-10 

Current (x 10 )A 


2.194 

3.194 

4.194 

5.194 

6.194 

6.294 

6.394 

6.494 

6.594 

6.694 

6.744 

6.794 

6.844 

6.894 
6.944 

6.994 
7.044 

7.094 
7.144 

7.194 
7.244 

7.294 
7.344 

7.394 
7.444 

7.494 
7.544 

7.594 
7.644 

7.694 

7.744 

7.794 

7.844 

7.894 

7.994 

8.094 

8.194 

8.694 

10.194 

11.194 

12.194 

13.194 
13.694 


4.030 

4.030 

3.990 
4.020 
4.010 
4.010 
4.010 

4.020 
3.980 
3.910 
3.810 
3.670 

3.480 
3,230 

3.020 
2.750 
2.620 
2.440 
2.285 
2.120 

1.990 
1.920 
1.855 
1.800 
1.660 
1.700 
1.690 
1.650 
1.640 
1.605 
1,600 
1.570 
1.560 
1.560 
1.540 
1.540 
1.505 
1.490 

1.490 

1.470 

1.480 
1.470 
1,450 


Concentrati on 
(std cm^/cm^) 

.0000 
.0000 
.0040 
,0020 
.0020 • 
.0020 
.0020 
.0010 
.0050 
.0121 
.0226 
.0379 
.0599 
.0914 
.1204 
.1618 
.1836 
.2162 
.2468 
.2824 
.3130 
.3305 
.3476 
.3626 
,4035 
.3914 
.3944 
.4066 
.4097 
.4208 
.4224 
.4322 
.4355 
.4355 
.4422 
.4422 
,4541 
.4594 

.4594 

.4665 

.4629 

.4665 

.4737 



Run B4 
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Vapor core temperature: 74.85° 

Sink temperature: 40.1° 

Heat input: 9.653 W 

Calculated gas content: 4.309 std cm^ 

Zero-absorption current: 4.03(x 10"^^) A 


Axial position (cm) Current (x 10~^^)A 


2.194 

4.030 

3.194 

4.030 

4.194 

3.980 

5.194 

4.010 

6.194 

3.990 ■ 

7.194 

4.020 

7.694 

4.020 

7.794 

4.000 

7.894 

3.990 

7.994 

3.990 

8.094 

3.980 

8.194 

3.870 

8.244 

3.820 

8.294 

3.630 

8.344 

3.400 

8.394 

3.140 

8.444 

2.870 

8.494 

2.750 

8.544 

2.490 

8.594 

2.295 

8.644 

2.165 

8.694 

2.020 

8.744 

1.845 

8.794 

1.740 

8.844 

1.160 

8.894 

1.560 

8.944 

1.150 

8.994 

1.460 

9.044 

1.425 

9.094 

1.405 

9.144 

1.395 

9.194 

1.370 

9.244 

1.360 

9.294 

1.350 

9.344 

1.320 

9.394 

1.320 

9.494 

1.310 

9.594 

1/310 

9.694 

1.290 ' 

10.194 

1.270 

10.694 

1.270 

11.194 

1.260 

12.194 

1.230 

13.194 

1.220 

13.694 

1.220 


Concentration 
(std cm^/cm^) 

.0000 

.0000 

.0050 

.0020 

.0040 

.0010 

.0010 

.0030 

.0040 

.0040 

.0050 

.0163 

.0216 

.0424 

.0697 

.1036 

.1428 

.1618 

.2068 

.2447 

.2723 

.3057 

.3503 

.3796 

.4192 

.4355 

.4524 

.4701 

.4829 

.4904 

.4942 

.5038 

.5078 

.5117 

.5238 

.5238 

.5279 

.5279 

.5362 

.5447 

.5447 

.5490 

.5622 

.5667 

.5667 



Run B5 
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Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


78.87°C 
40.22°C 
11.916 W 
4.195 std cm^ 

1 f) 

4.03{x 10" ) A 


Axial position (cm) 


Current (x 10~^^)A 


Concentration 
(std cm^/cm^) 


2.194 

4.020 

.0010 

3.194 

4.030 

.0000 

4.194 • 

3,980 

.0050 

5.194 

4.020 

.0010 

6.194 

4.010 

.0020 

7.194 

4.010 

.0020 

8.194 

3.990 

.0040 

9.194 

4.020 

.0010 

9.294 

3.970 

.0060 

9.394 

3,830 

.0205 

9.444 

3.670 

.0379 

9.494 

3,380 

.0722 

9.544 

3.120 

.1063 

9.594 

2.810 

.1522 

9.644 

2.540 

.1977 

9.694 

2.330 

.2376 

•9.744 

2.120 

.2824 

9.794 

1.935 

.3267 

9.844 

1.790 

.3654 

9.894 

1.673 

.3995 

9.944 

1.560 

.4354 

9.994 

1.440 

.4773 

10.044 

1.350 

.5117 

10.094 

1.276 

.5421 

10.144 

1.230 

.5622 

10.194 

1.180 

.5850 

10.244 

1.175 

.5874 

10.294 

1.166 

.5916 

10.344 

1.140 

.6042 

10.394 

1.130 

.6091 

10.444 

1.125 

.6116 

10.494 

1.120 

.6141 

10.594 

1.110 

.6191 

10.694 

1.080 

.6345 

10.794 

1.080 

.6345 

10.894 

1.080 

.6345 

11.194 

1.070 

.6398 

12.194 

1.080 

.6345 

13.194 

1.135 

.6066 

13.694 

1.030 

.6607 



Run B6 


^3 


f 

I 


Vapor core temperature; 83.1°C 

Sink temperature: 40.25°C 

Heat input: 14.336 W 

3 

Calculated gas content: 4.157 std cm 

”10 

Zero-absorption current; 4.04(x 10 )A 

Concentration 
3 3 

Axial position (cm) Current (x 10 ^>)A- (std cm /cm ) 


2.194 

4.040 

.0000 

3.194 

4.030 

.0010 

4.194 

4.030 

.0010 

5.194 

4.030 

.0010 

6.194 

4.020 

.0020 

7.194 

3.980 

.0060 

8.194 

4.000 

.0040 

9.194 

4.040 

.0000 

9.694 

4.030 

.0010 

9.994 

3.990 

.0050 

10.094 

3.980 

.0060 

10.194 

3.980 

.0060 

10.294 

3.890 

.0152 

10.344 

3.830 

.0215 

10.394 

3.690 

.0367 

10.444 

3.470 

.0622 

10.494 

3.050 

.1172 

10.544 

2.780 

.1580 

10.594 

2.420 

.2211 

10.644 

2.180 

.2702 

10.694 

1.970 

.3191 

10.744 

1.790 

.3666 

10.794 

1.605 

.4220 

10.844 

1.490 

.4607 

10.894 

1.370 

.5051 

10.944 

1.270 

.5460 

10.994 

1.165 

.5935 

11.044 

1.080 

.6359 

11.094 

1.030 

.6620 

11.140 

0.980 

.6673 

11.194 

0.965 

.6990 

11.244 

0.936‘ 

.7168 

11.294 

0.938 

.7155 

11.394 

0.908-' 

.7348 

11.494 

0.895 ■ 

.7434 

11.594 

0.890 ' 

.7468 

11.694 

0.890' 

.7468 

11.794 

0.880 

.7536 

12.194 

0.880>' 

.7536 

12.694 

0.880'. 

.7536 

13.194 

0.880 

.7536 

13.694 

0.880 

.7536 
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Run B7 


Vapor core temperature: 91.3^C 

Sink temperature: 40.2°C 

Heat input: 18.807 W 

3 

Calculated gas content: 3.915 std cm 

-10 

Zero-absorption current: 4.03(x 10 ) A „ ^ j.. 

Concentration 

Axial position (cm) Current (x 10 ^?-)A (std cm /cm ) 


2.194 « 

4.020 

.0010 

3.194 

4.020 

.0010 

4.194 

2.980 

.1263 

5.194 

4.020 

.0010 

6.194 

3.980 

.0050 

7.194 

3.980 

.0050 

8.194 

3.940 

.0091 

9.194 

4.070 

-.0039 

10.194 

4.050 

-.0020 

10.694 

4.030 

.0000 

11.194 

4.030 

.0000 

11.294 

4.000 

.0030 

11.394 

3.900 

.0132 

11.494 

3.900 

.0132 

11.594 

3.880 

.0152 

11.694 

3.870 

.0163 

11.794 

3.640 

.0413 

11.894 

3.300 

.0823 

11.994 

2.680 

.1733 

12.094 

2.180 

.2690 

12.194 

1.460 

.4701 

12.294 

1.180 

.5850 

12.394 

0.985 

.6859 

12.494 

0.860 

.7661 

12.594 

2.741 

.8604 

12.694 

0.650 

.9487 

12.794 

0.650 

.9487 

12.894 

0.640 

.9594 

13.194 

0.630 

.9704 

13.694 

2.610 

.9932 
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APPENDIX B 


SELECTED TABULATIONS OF CALCULATED TEMPERATURE AND 
CONTENTRATION PROFILES BASED ON THE 1-DIMENSIONAL 
INSOLUBLE GAS MODEL 


NOTE: The sign appearing in the printout is a misprint for the;numeral "3". 



Run A1 


66 


Vapor Temperature 49.65°C 

Sink Temperature 14.0°F 

Gas Content 4.342 std cm^ 

Calculated Heat Transfer Rate 20.47 W 


Axial Position 
(cm) 

T"^ 
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Run A1 (continued) 



.261957 

. .26669 

. 23 33 

? . 1 r-4 1 ? 

.24^283 

. 240404 

. 24 34 35- 

2.^'05t.l 

. 227357 

. 2- 1285 

. 247245 

?:..?3ri 

..212347 

. 21 5267 

. 25 0751 

’ . 56553 

, 1 3 ,■’9 0 1 

. 200296 

. 25 33.34 

i. 3'"tOC"r. 

. 124172 

. 186=83 

. 252966 

£. 33154 

. 171411 

, 1 73193 

. 2597 1 5 

£'. 36 5 03 

. 15947 

. 160379 

. 262257 

■5.3 "*451 

. 14.3305 

. 149572 

. 2646 01 

-■.4 55-=<9 

. 157S71 

. 1 ?992’=^ 

■ c ‘Zf *Z‘ 1 to 

-.45747 

. 123125 

. 1230U1 

. 26-S767 

. 43'5 36 

. 1190=7 

. 1 13745 

.270615 

5. 53044 

. 11053? 

.111119 

..='72324 

3. 55195 • 

. 1 0262 

. 1G30E2 

■ .=75303 

3.5 5-34 

. 95=*53?5- 1 

.?55=»72E-1 

. 275363 

-3. 6 14.? 3 

.833605E-1 

. 336297E-1 

.276714 

?■ 64637 

.fil«1536f-l 

.221457E-1 

. 277962 

?. 67785 

.759981 E-1 

.7611336-1 

. 279117 

■'0=‘55 

. 704357E-1 

. 7050466-1 

. 28 0186 

3. 74 052 

. 6526.99E— 1 

.b52?01£-l 

.281174 

- • "Fes 

.604654E-1 

. 604443E-1 

. 232088 

3. c 0372 

. 5.599S6E- 1 

.559427E-1 

.282934 

? . $ 5:52^ 

,51 S47 riE" 1 

.517622E-1 

. 28 37 16 

3. 3t.^-75 

,4=’=‘?C5E“1 

.473312E-1 

. 294439 

3. 39?c3 

. 444085E-1 

.44.2794H-1 

.2951 08 

3.93971 

. 041033 

, ijii4 0 ? 

, 285727 

•?. ^iA 19 

. 3 79966E— 1 

. 37.=:?9ie-l 

. 286299 

?. 99? 6.3 

. 351 331E-1 

. 0*34966 

. 986? 28 

5‘. 024 16 

.'324772E-1 

. 39-021E-1 

.287317 

3. 0555.4 

.30 t.t 1 -47 E— 1 

. 299359E-1 

. 237769 

3. 113713 

.277923E-1 

.27550=6-1 

. 2.5812= 

11261 

.256174E-1 

. 254 6496 - i 

. £‘ 5 5 ^ r f 

;*. 15 0 06 

. 2325841=- 1 

. ■ 2.4 I't ;.t — 3 

. c 3 r'3 r'9 

; . 1 3 3 5 r 

. 21 3449 .--: 

. 23 6 -.■-.26-'. 


:'.,2i3 06 

. 201653E-1 

. 1 'r 2 27 4- — 1 

. 2? 

:-4454 

. IS' 1 1 :-2-i 

. U ’'=427 

m cl '1* y 

1'. r’i’3.1-‘c' 

. 1 7 1 7 i 

. 1 7 0 ‘J b.— 1 

.23 01 U=' 

3. 3075 

. 0 1 5044 

. l5r-,780 6-l 

. 23 0?4if 

3. JiV-i 

. 1 4t-l 4?t - 1 

. 144-4 3-=-1 

. 05- r 

3. :-.7P4r 

. 1 34794=- 1 

, : 3 5 ?4 26- 1 

.2307 67 

3.4f|l-^'^ 

. 1.24295S-1 

. 12260= 6-1 

. 2 =‘0955 

J • -4 :;I r* 

. 1 1 4536b- I 

. 1 1 :nb46-i 

. 2‘3U2? 

•1 . 46.^97 

. 1 056 = '"'?— 1 

. ui 04-= r 

. =3125-' 

. 4 32.4 

. ?r?.5’37b“ 2 

. 96 06«E- 2 

. 2'31-^ 36 

J« “'■? 

. 227259E— = 

. .?.?4 3 3.iE-2 

. 5 315 71? 

3 . 559 36 

. 3=2 5?b-? 

, 915181 £-2 

. 29 1699 

- 1 . 59 '155 

.7t.197.lb-2 

.75070.19-2 

. 2’=' 1 9 1 

-1 . 2-c ? ?3 

. 7 019466:- 2 

.6914726-2 

, 29 1 9 1 =^ 

V 


ORIGINAIi PAGE IS 
05} POOR QUALHY 


68 


Run A1 (continued) 
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Run A3 


Vapor Temperature 
Sink Temperature 
Gas Content 

Calculated Heat Transfer Rate 


Axial Position 
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Run A3 (continued) 
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Run B1 


Vapor Temperature 63,35°C 

Sink Temperature 40.12°C 

Gas Content 4.49 std cm^ 
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Run B1 (continued) 
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Run B1 (continued) 
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Run B2 
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Gas Concentration 
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Run B2 (continued) 
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Run B2 (continued) 
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.9471 j lE-3 

. 34996 1 E-9 

• 39774^1 

7 . 1) " '?4 

. .599 1 75E- 3 

. 395944E-9 

. 397797 

1 ■ 3 4 4 1 

. 31 9^00 E- 9 

. 90903 4E-5 

. 397844 

7.665 

. 19 79 3 PE- 3 

, 3 6^33 

. 397887 

7. 665'5? 

• 1 849 3 

. 130897E-3 

. 3979 36 

7. 76605 

.1718 ?’PE- 5 

. 168795E-P 

. 3^7969 

7. 756 57 

. 16071 8E-3 

. 157717E--9 

. 397994 

7.79711 

. 150648E-9 

. 14779.5E— 9 

. 398093 

7. 31 76*^ 

. 141*:>149 — 

. 1 33903E-P 

. 39c‘049 

7 • 64617 

. 135556E-9 

. 1 if 097 IE -9 

. 39807= 

7 . 97'5‘7 

. 196419E-9 

« 19^4 4 ^ ^ 

. 59'= 099 

7. ?09.?3 

. 130157E-9 

. 11778 9E-3 

. P?811 

7. 959 75 

. 114796E-9 

. 119449E-3 

. 39E!l £t 

7 . ‘P7 069 

. 11 0093 E- 9 

. 1 0789 3E-9 

. 398139 

3. 00091 

. 1 0699:-£-9 

. 1 0409E-9 

. 39515 

3. 031 34 

. 103099E-9 

.101 013E-9 

. if9Sl 59 

5. 06 197 

. 1 00673E-3 

. 9864 1 7E— 5 

. 3931 66 

6 • 0994 

, 9S9c:.6 3t- 3 

. 64.-=iS>P4P--' 

. 39? 171 

?. 19 39? 

. 9794 36 E- 3 

# r* ^*4 *r*H“ n* 

. 393174 

5. 15346 

. 976034E-3 

. 956205e-3 

. 398175 
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APPENDIX C 

LISTING OF COMPUTER CODE FOR 1-DIMENSIONAL 
MODEL OF CONDENSATION HEAT TRANSFER IN THE 
PRESENCE OF AN INSOLUBLE GAS 


NOTE: The sign appearing in the printout is a misprint for the numeral "3". 



! T 79 

- ’SI Oo-FEB-76 

- L s 1 ■ ' ■ > , y i;£g4i , C <2g4> > D ^a44> . P i'a64> 

. -■ .-hSl fT)=fl4*t/!P(OS>*T»*'=)3 

' ’ FNPn<T'«=qi*ESP<'4c-'T'' 

: =il) ZO SraTPt CONDEMiEP LENSTH 

- "-iTl=l 16. 

• :»=ID SI. -3£* nil j Ilia H-H8U_ 

'I '.8T8 .41 ?6. . 0??5?. . 01 1?- -4?l>? 

READ ■n.Va.*'3. !fl-VfiPOR< CONST. IN C»D,riT. OF Wfipof;. 

=E tiOTfl 1-047.P.31E-3. 51 000. 

REflri 80,85,82 !LN’Pl= Rij + + 8s'*lN^T5 <P IN RTMOSPHEPES) 

16 tiflT« ?0.5S75,-4344.47>-2.?0124 

40 PEflD N-55-M !M=SRS-Z0HE IMCF . . Za=DIM. LENGTH. M=INCR. IM FREE Z. 

41 DftTR 114., e. 5. 10. 

45 PER E5 >MRX. EPPO® PERMITTED IN T-',I8lL PER ITERATION 

46 D8TR , OOl 

51 INPUT"T-nrpgp=- ;t1 

52 INPUT •T-SIflK='"T2 

5;- INPUT STD. CC. OF GRS= SG^ 

6" H=Za-'N:tU=rt+l:NE=N-n5; ril=s1+l : "lE=ri+N: M3=M+N+1 
o5 81=57.5. 15^E!>'P I flu'' s pi==NFii-'T 1+57?. 15., :fl4=Rl^Pl 

70 Um=T 3+E73. 15:U3=Tl-n:C5=LI?/'4. : C 2=H**£*Gl*ijH»li‘5'' GE^Vl»''i*VlJ> 

71 ri = £PP<IVU*M£>'S£., : G?=l . Cl . .-'51 + 1 . '■ ''13=63^55 

7^ FI =— c. — Cl , +>51 /Sa> E! t5=“Gl E., G 2 

7? H0=Z1+H 
90 MAT lil=';ON<'in 
^1 MAT M= *1,1 

=5 'M’=.P95:>: >,1>=0. ; pa,="> n IROJUIT for calc, of C'Z 1 ,Di.£) 
luo FOP 1=3 TO 

105 7> I'i=£VP>'-i'I-l Ji+rij !M I I+n)=i,i5*'wl' 5NE1‘‘T I 
nil pt£i=i./'Di 

115 FOP 1=2 TO M3 iP-riRTRiy FOR 'JflLL TEMP. SOLUTION 
120 P<l)=l. /-CPi-F . I-D I : NEsT I 
15u T3=U. : GD TO ZGO 
130 T5=0. 

200 FGP I=? to Ml JCRLC. C AMD D URTPIV? FOP I'RPOr TEMP. SOLN. 

205 T=Un+ua+v>. I) ;Ct=fl2'T:''=fl4*E'CP<'CE i+T^+Rs 1C7=R3-CG 
21 0 ?=C5+'"iI + l >-V> I-m ♦<'37.' a. -7->+.;6/C7-,5 j ^T 
215 F=C2+?0P>'T1 + C1.'''‘.-1. ’.'C7 
=?n b=-i'e.+f;i_^i._E''+c ci-n 

235 Ccl ■ = Cl.+E' .-E: DCI< =->.F+i.i>'I+.’ ■ + Cl . -E ■ *D <I - 1 • • 'B 
240 MEv-T I 

=50 VCMl ■ = >'F+i,H'ri?'>+2. ♦DCN,'' <2.+F+E.+CCN> ' 

2t;0 FGP I=N TO ? STEF -I • -'pprip TEflPEPRT'JPE PROFILE 
2^5 •‘i'n=P(I>-C ‘'P*“<1 + 1'' : NEXT I 
230 D 2> = CE2-i>I3'7E1 

255 FOP 1=2 TO M !D MRTPI'; FOP i.ILL TEnpeFATUPE CALCS. 

290 DCI' = CP2-DcI-ij)j *p<n :NEVT I 

s uo FOP i=ni TO n= 

?05 D<I'' = CE2*>'sI“m-DCl-l'"+P<'D iNE^iT I 

210 M>M3'' = 'B£+'J'UU.>-2.+D'M2'i J/'iIl-S.+F <i>121 J :M4=M2 

325 FOR I=M2 TO 2 STEP-I 

330 T1='.|CI% jl,lcn=0Cl)-P IT+HCI + n ;T2=ATSCT1-i>)CI’> 

S‘35 IF T2<T? GO TO 34 : 

'40 T3=T2 
545 s(FPT I 

400 IF T?'E5 GO TO 1=0 

600 C3=cpi-FNFOCiJ0)> -'UO* T=Uii+LIP+'' ■'!> s C cn = CPl-FMPO CT, ’''T 
605 T='JO+U9 *''Cni>!C*'N 1'=>,P1-FNP0CT^ > ''T 
610 5=.5+>'c<n+C ('fu,'>: C'2=.5+> MCI ) +i,i>'ri2i ' 

612 FOR 1=2 TO ’ll ! IMTE'EFATE rlCG AND TEMP. PROFILES 

613 r(2=C'2+I.ICr) : N£>;t I 

615 FOE 1=2 TO M ‘INTEGRATE GAS AND TEMP. PROFILES 

6-0 T=UO+U?*''ct'i ; C cI'=>'pi-FNPO>'T'>i ■'T 

625 G=>3+C I): C‘2=‘''2+M Cl+t)'' : NEvy I 

o30 6=V1+H0*G: n2=U3+G2+HO*02: Z7= ‘'GO-G) f''l+C 2 J 

1 4 0 Z5=Z 0- ''Z7+Z2+Z 1 +H u ♦M'> : r> 1 =U =+GE+U < 1 '> +Z5 : A 2 =■', 1 -^02 

645 fqp 1=1 TO M3 

650 D'I)=Z5+H0* I-l •NE>{T I 

670 PPIMTS PRINT- CDN'SENTPATIQN AND TEMPERATUPE PROF ILE3“ ! PRINT 

671 PRINTspciNT Ml . ^12 . A3: PF INT 
620 FQP 1=1 TO M 
6^0 FcinT DcD .iJCIj iNENT I 
■?^1 fqp j=fit TD ft? 

692 FPINT D.I) .MsD ,V- I-Ml .C> I-Mn NE-^T I 
70U F'3P 1=1 TO N1 
710 D<I'i=HO+CI-l):MEPT I 
720 PRINT; PRINT 
72f FOP I=l TO N1 

740 PPIMT CIJ ,MiM+I> , >/CD ,CCl-> :NE^T I 
300 MNT ppIMT 
3110 FND 


OElGINAi; 

OF POOR QUA3jin3 



